Persistent superfluid phase in a three-dimensional quantum XY model with ring 

exchange 
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We present quantum Monte Carlo simulation results on a quantum S — 1/2 XY model with ring 
exchange (the J-K model) on a three-dimensional simple cubic lattice. We first characterize the 
ground state properties of the pure XY model, obtaining estimations for the energy, spin stiffness 
and spin susceptibility at T = in the superfluid phase. With the ring exchange, we then present 
simulation data on small lattices which suggests that the superfluid phase persists to very large 
values of the ring exchange K, without signatures of a phase transition. We comment on the 
consequences of this result for the search for various exotic phases in three dimensions. 



The search for microscopic spin 1/2 Hamiltonians that 
exhibit a spin liquid phase continues with the hope that 
these studies will provide insight into the type of mate- 
rials that can support such a phase. Focus has recently 
turned to multi-particle ring exchange models^^ which 
have been identified as promising candidates for the pres- 
ence of spin liquid as well as other exotic phases and 
critical points. One such Hamiltonian is an easy-plane 
5=1/2 J-K model on the square lattice, with an XY 
exchange coupling J and a four-site ring exchange cou- 
pling K. Quantum Monte Carlo studies 1 - of this model 
have shown that, at a critical value (K/J) c , there ap- 
pears to be a continuous quantum critical point sepa- 
rating an XY-superfluid phase from a valence bond solid 
(VBS) phase. Recently, Senthil et alm> have proposed that 
at this type of critical point - separating conventional 
phases described in terms of different broken symmetries 
- one can have an isolated spin-liquid point character- 
ized by an emergent global U(X) symmetry and decon- 
fined spinon excitations. In two dimensions, the theory 
finds that the U(l) spin liquid is unstable away from the 
critical point because of the proliferation of monopolc 
excitations, which are only absent precisely at (K/J) c . 
They suggest that it is this instability that is respon- 
sible for the formation of the VBS phase seen in the 
simulations 1 for K/J > (K/J) c . However, in three di- 
mensions, monopole excitations can be gapped, so within 
the theoretical framework^ an extended region of stable 
U(l) spin liquid phase could exist £ Motivated by this 
possibility, we have carried out a quantum Monte Carlo 
study of the simple cubic J-K model to look for signa- 
tures of a zero-temperature phase transition out of the 
superfluid phase (at large J/K) to a quantum disordered 
or insulating phase (at large K/J). 

The Hamiltonian under study is 



H=-J^B ij -K £ P ijkl , 

(ij) (ijkl) 



(1) 



where bond and plaquette exchange operators are defined 
as 



Here, (ij) denotes a pair of nearest-neighbor sites and 
(ijkl) are sites on the corners of a plaquette (see Fig.^l. 
In three dimensions, each spin is shared by six bonds 
and twelve plaquettes, which form the faces of the simple 
cubic lattice. For K = this is the standard quantum 
spin- 1/2 XY-model, or, equivalently, a hard-core boson 
model at half-filling, studied previously using quantum 
Monte Carlo techniques by Pedersen and Schneider £ The 
model is analogous to a quantum lattice gas model of 
He 4 , and undergoes a lambda-transition at T/J w 2 to 
a bulk superfluid statei£*£ The if-term corresponds to 
retaining only the x- and y-spin operators of the full four- 
spin cyclic exchange studied in relation to some high-T c 
parent compounds. 8 
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FIG. 1: The labeling of the plaquette indices in Eq. Q are 
illustrated on left. At right, three intersecting plane segments 
of the cubic lattice show that each spin site shares six bonds 
and twelve plaquettes with other sites. 

Our simulations use the stochastic series expansion 
(SSE) quantum Monte Carlo method)2ii2, which is based 
on importance sampling of elements of the power se- 
ries expansion of the partition function. The method 
is numerically exact within statistical error bars, and has 
broad applicability to a wide range of spin and boson 
models where the sign problem can be avoided. In order 
to efficiently sample off-diagonal J terms in the Hamil- 
tonian, the recently developed "directed-loop"^ 1 * 1 ^ algo- 
rithm is used. However, for large values of the param- 
eter K, the loop moves become ineffective, and a new 
sampling procedure was developed^ This new "multi- 
branch cluster" update is a Swendson-Wang or Wolff- 
type sampling scheme generalized to flip clusters of pla- 
quette operators in the SSE operator- list. It is observed 
to significantly decrease autocorrelation times of the sim- 
ulations for large K, even for K/J — > oo in two dimen- 
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sions. Details of the SSE algorithm, simulation procedure 
and calculation of physical observables can be found in 
Ref . E3 

The energy per spin is simply calculated in terms of the 
SSE estimator for the order of the power series expansion, 
n, and the temperature: 



E 



(n)T 
N ' 



(4) 



where N = L 3 is the system size. The main quantity of 
interest in this paper is the spin stiffness, 



d 2 E((/>) 
d<t> 2 



(w 2 ) + (w 2 )) 



3N 



(5) 
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where is a twist imposed on bonds in either the x, y or 
z lattice directions. In Eq. © we calculate the deriva- 
tive at <p = using the winding number estimators W XjVtZ 
which measure the net spin current across the periodic 
boundaries in each directioniiflii^ We have also calculated 
in some cases the spin susceptibility x, easily obtainable 
in the S z basis using the squared magnetization estima- 
tor, 



- TN (\Y. S 



(7) 



In three dimensions, even the smallest useful lattice 
sizes have Hilbert spaces far too large to allow an exact 
diagonalization study to be performed. In order to test 
the Monte Carlo and to characterize the K = ground 
state, we study the asymptotics of the pure XY model, 
setting J = 0.5 in the Hamiltonian Eq. QJ. Second- 
order spin-wave calculations carried out by Weihong et 
alH give numerical values for the ground state energy 
per spin and spin susceptibility as Eq w —0.7914 and \ ~ 
0.1517, without error bars. Although we are not aware of 
any detailed numerical results for the zero-temperature 
spin stiffness, mean field theory predicts a value of p s = 
S 2 = 0.25 using Eq. JSJ. One might expect that quantum 
fluctuations slightly increase the value of p s in the full 
5 = 1/2 XY model, as observed in two dimensions^ 

We simulated the pure XY model (K = Q, J = 0.5) for 
several lattice sizes, L= 6, 8, 10, 12, and 14, using the full 
SSE code, for several values of T/J down to 0.04. The 
relevant estimators either showed absolute temperature 
convergence to within error bars, or were extrapolated 
to T = 0. Data obtained from the SSE simulation for 
the ground-state energy per spin, Eq, is most accurate. 
Chiral perturbation theory^ predicts that the energy de- 
pends on system size and spatial dimension (d) as 



E (L)-E 



C 
Ld+T> 



(8) 



where C is some constant. This finite-size dependence 
has been shown to be in good agreement with similar 
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FIG. 2: Ground-state energy versus 1/L 4 for system sizes 
L = 6,8,10,12 and 14. The curve is a fit to Eq. ©. 



quantum Monte Carlo studies of the two-dimensional 
(2D) XY models Setting d = 3 in Eq. ©, we use 
our data to obtain an extrapolation to L — ► oo, giving 
E = -0.79182(2). The fit is illustrated in Fig.|3 which 
has a minimum in the chi-squared value per degree of 
freedom for C = —1.14. Data for the spin susceptibil- 
ity and spin stiffness have considerably larger error bars, 
within which the larger system size data is converged, 
making a finite-size scaling extrapolation irrelevant. Our 
best estimates for the ground state values in these cases 
are x = 0.1500(2) and p s = 0.2623(2) 
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FIG. 3: The finite-temperature behavior of the spin stiffness. 
Data for both L = 6 and L — 24 is included for K/J = to 
illustrate finite-size effects at the lambda-transition. 

Turning to the full model, we performed simulations 
including the K term in Eq. QJ, keeping J fixed at 1/2. 
First, the T — > temperature convergence was estab- 
lished for each value of K/J simulated. Fig. |21 compares 
the finite-temperature behavior of the spin stiffness for 
an L = 6 system at K/J = 20 versus the XY-point at 
K/J = 0. Results for the ground-state spin stiffness ver- 
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sus K/J are illustrated in Fig. 21 Data for L — 4 and 
L = 6 up to K/J = 48 was taken at temperatures below 
T/J = 0.5, depending on the specific K/J where good 
T — > convergence was established. However, the three 
largest K/J values for L = 6 began to show sticking, 
evidence that the Monte Carlo time-scales are getting 
very large or ergodicity is being lost. To combat this, 
data for L — 8 was taken at slightly higher temperature, 
T/J < 2.0 (for large K/J), which was nonetheless within 
the range of the zero-temperature convergence. 

For the 2D J-K model, the zero-temperature critical 
point separating the superfluid phase from the VBS oc- 
curs at (K/J) c « 8,i However, as illustrated in Fig. ^ 
for the three-dimensional (3D) case the superfluid phase 
is clearly much more robust. In particular, extensive 
simulations for the L = 6 system, out to K/J = 144 
fail to show an abrupt transition out of the superfluid 
phase. Similarly, data for L = 8 remains strongly su- 
perfluid for K/J up to at least 60. As evident from 
Fig. finite-size effects for the spin stiffness become 
much more pronounced as K/ J gets large. Further, ob- 
servation of the staggered magnetization showed no indi- 
cations of the onset of any Neel (or boson charge-density 
wave) ordered state at any K/J which was simulated. 
Unfortunately, algorithm performance becomes less effi- 
cient for this model for large values of K/ J. This, cou- 
pled with the fact that the lattice-plaquette coordination 
contributes an extra factor of three in the CPU time size 
scaling of the algorithm in three dimensions as compared 
to two dimension, makes a more detailed study of larger 
lattice sizes impractical at this time. 
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FIG. 4: Ground-state spin stiffness for the 3D J-K model. 
The XY exchange coupling was set to a constant, J = 0.5. 
Open symbols indicate data which showed signs of ergodicity 
loss, and therefore error bars may be under-estimated. 



The consistency of data for the two largest system 
sizes points toward robust superfluidity in the system well 
above the value for (K/J) c where the onset of insulating 
behavior occurs in the 2D modeli Although impossible 
to prove without a detailed finite-size scaling study, we 
believe that the data supports the conjecture that the 
superfluid ground state persists in the model, becoming 
critical only at K/J —> oo. This of course leaves the in- 
teresting question regarding the nature of the K/ J — > oo 
phase open for future study. One suggestion is that the 
J = model actually has long-range in-plane ferromag- 
netic order, in accord with a simple spin-wave analysis. 17 
The critical value for the spin stiffness in this case would 
be p s = 0, however for finite J as in Fig. 0] p s may be 
non-zero in the thermodynamic limit. 

In summary, we have presented data from SSE quan- 
tum Monte Carlo simulations of the 3D simple cubic J-K 
model. Simulations of the XY point reveal a ground state 
energy, spin susceptibility, and spin stiffness of 

E = -0.79182(2), 
X = 0.1500(2), 
p s = 0.2623(2), 

consistent with predictions. In particular, the ground 
state energy agrees with second-order spin- wave theorj*i4 
to within 0.05%. The full J-K model appears to remain 
superfluid for all finite K/ J, with no phase transition 
signature up to K/J = 144. The large K insulating 
phase is still unknown, however no indications of cither 
a charge-density wave or any quantum disordered state 
are present. The robustness of the superfluid phase is 
consistent with the notion that the relative strength of 
quantum fluctuations introduced by the ring-exchange 
(K) term is reduced as one increases the dimensionality of 
the system. Finally, the failure of this model to support 
an emergent U(l) spin liquid illustrates the difficulty in 
realizing exotic phases in microscopic quantum models 
suggested by arguments arising from generalized actions. 

Acknowledgments 

The authors would like to thank T. Senthil for moti- 
vating discussions. We appreciate extensive and help- 
ful collaboration with A. Sandvik. Thanks also to 
L. Balents, A. Burkov, M. P. A. Fisher, M. Hermele, 
and O. Motrunich for useful conversations. Supercom- 
puter time was provided by NCSA under grant number 
DMR020029N. This work was also supported by the De- 
partment of Energy, Grant No. DOE85-45I97. 



1 A. Sandvik, S. Daul, R.R.P. Singh, and D.J. Scalapino, 
Phys. Rev. Lett. 89, 247201 (2002). 

2 L. Balents, M.P.A. Fisher, and S.M. Girvin, Phys. Rev. 



B 65, 224412 (2002); V. Rousseau, G. G. Batrouni, 
and R. T. Scalettar, cond-mat/0312018 M. Hermele, 
M .P. A. Fisher, and Leon Balents, Phys. Rev. B 69, 



4 



064404 (2004); V. Gritsev, B. Normand, and D. Baeriswyl, 
Phys. Rev. B 69, 094431 (2004). 

3 A. Paramekanti, L. Balents and M. P. A. Fisher, Phys. 
Rev. B 66, 054526 (2002). 

4 T. Senthil, A. Vishwanath, L. Balents, S. Sachdev, 
and M. P. A. Fisher, Science 303, 1490 (2004); 
T. Senthil, L. Balents, S. Sachdev, A. Vishwanath, and 
M. P. A. Fisher; cond-mat/0312617 

5 D. A. Huse, W. Krauth, R. Moessner, and S. L. Sondhi, 
Phys. Rev. Lett. 91, 167004 (2003); R. Moessner 
and S. L. Sondhi, Phys. Rev. B 68, 184512 (2003); 
O. I. Motrunich and A. Vishwanath, cond-mat/0311222 
O. I. Motrunich and T. Senthil, cond-mat/0407368 

6 M. H. Pedersen and T. Schneider, Phys. Rev. B 53, 5826 
(1996). 

7 D. D. Betts, in Phase Transitions and Critical Phenomena, 
edited by C. Domb and M. S. Green (Academic, London, 
1974), Vol. 3, pp. 569-652. 

8 M. Roger and J. M. Delrieu, Phys. Rev. B 39, 2299 (1989); 
Y. Honda, Y. Kuramoto, and T. Watanabe, Phys. Rev. B. 



47, 11329 (1993). J. Lorenzana, J. Eroles, and S. Sorella, 
Phys. Rev. Lett. 83, 5122 (1999); R. Coldea, et al. Phys. 
Rev. Lett. 86 5377 (2001); B. Normand, and A. M. Oles, 
cond-mat/0403481 
9 A.W. Sandvik and J. Kurkijarvi, Phys. Rev. B 43, 5950 
(1991); A.W. Sandvik, J. Phys. A 25, 3667 (1992). 

10 A.W. Sandvik, Phys. Rev. B 56, 11678 (1997). 

11 A.W. Sandvik, Phys. Rev. B 59, R14157 (1999). 

12 O.F. Syljuasen and A.W. Sandvik, Phys. Rev. E 66, 
046701 (2002). 

13 R. G. Melko and A. W. Sandvik, in preparation. 

14 Z. Weihong, J. Oitmaa, and C. J. Hamer, Phys. Rev. B 
44, 11869 (1991). 

15 A. W. Sandvik and C. J. Hamer, Phys. Rev. B 60, 6588 
(1999). 

16 P. Hasenfratz and F. Niedermayer, Z. Phys. B 92, 91 
(1993). 

17 L. Balents, private communication. 



